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Abstract 

It is shown that linearized gravitational radiation confined in a cavity 
can achieve thermal equilibrium if the mean density of the radiation 
and the size of the cavity satisfy certain constraints. 

Is it possible, in principle, to confine gravitational radiation in a box, for 
a time long enough for it to achieve thermal equilibrium? This issue has been 
investigated before by various authors [1, 2, 3], not all of whom have arrived 
at the same conclusion. The purpose of this note is to pursue a particular 
line of thought which does not seem to have been looked at in earlier work. 
We argue that it is possible to thermalize gravitational radiation, starting 
from a non-equilibrium configuration, so long as certain constraints (to be 
discussed below) are satisfied. 

Let us begin by summarizing the arguments given by Smolin [1]. He 
showed, by considering various mechanisms of absorption of radiation by 
matter, that no realistic material is an efficient absorber of gravitational ra- 
diation, except over a narrow band of frequencies. These mechanisms include 
absorption of classical gravitational radiation by classical matter, ionization 
of atoms by gravitons, and phonon excitation. He concluded that a state of 
thermal equilibrium between radiation and matter cannot be reached in a 
finite time. 
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Garfinkle and Wald [2] presented a counter-example to the result by us- 
ing a shell of charged matter balanced just outside its Schwarzschild radius. 
However Dell [3] later argued that such a shell is unstable in the presence of 
electromagnetic radiation. 

Whatever the nature of the confining cavity, the absorption coefficient / 
will always be smaller than one. If the box is made thick enough to raise / 
to one, it undergoes gravitational collapse. 

However, it is not absolutely essential, for the purpose of attaining equi- 
librium, to have / close to unity. For any cavity, it is possible to define a 
leakage time-scale ti over which, say, a fraction 1/e of the initial radiation 
escapes. Also, for an initial configuration p(y) - the energy-density of gravi- 
tational radiation in the cavity - an equilibrium time-scale t# can be defined, 
over which the radiation may achieve equilibrium. What needs to be settled 
is whether can be smaller than t^, for a general initial configuration p(u). 
We argue below that this is possible, and ts < need not imply gravi- 
tational collapse of the cavity or of the radiation inside it. Moreover, we 
have no need to concern ourselves with any specific mechanism for confining 
radiation. 

Our argument is analogous to the standard discussion of Brownian mo- 
tion of a molecule in a fluid, where one shows that as a consequence of the 
random collisions and viscous drag, the molecule attains a Maxwellian veloc- 
ity distribution. A variant of this argument was used by Einstein [4] to derive 
the Planckian distribution for electromagnetic radiation in equilibrium in a 
cavity. The argument went as follows. He assumed that the radiation was 
in equilibrium with a molecule (say, with two energy levels). The molecule 
absorbs and re-emits at a characteristic frequency uq. The momentum trans- 
fers to the molecule take place through random collisions with the quanta 
of radiation, and through the systematic drag force it experiences because 
it sees a Doppler-shifted radiation while in motion. Einstein showed that 
if the molecule has a Maxwellian velocity distribution, then the exchange of 
momenta with radiation preserves this distribution only if the energy density 
of the radiation obeys the Planck law. 

In his paper, Einstein concerned himself with the equilibrium situation. 
The discussion of momentum transfers to the molecule assumes greater sig- 
nificance in a study of approach to equilibrium. Now, there is a systematic 
evolution of the r.m.s. velocity of the molecule, which may be approximately 
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described by the Langevin equation 



M^ = -av + /3(f), (1) 

a being the drag coefficient and (5 the fluctuating force (a is positive). Stan- 
dard discussions [5] of Brownian motion show that the molecule attains a 
Maxwellian distribution over a time-scale 

M 

t E w — . (2) 

a 

In particular, the 'fluid' pushing the molecule around could be electro- 
magnetic radiation, and (2) will continue to hold. [The assumptions which 
go into using the Langevin equation as well as subtleties related to the de- 
pendence of the drag coefficient a deserve discussion; this is presented in the 
last section since they are unrelated to the main theme of this paper.] Once 
the Maxwellian distribution has been reached for the molecule, Einstein's 
argument can be applied to arrive at Planck's law for the radiation. 

We now return to the case of linearized gravitational radiation in a cavity, 
interacting with a two-level atom. 

To begin with, we compare this system with the one consisting of an atom 
in electromagnetic radiation. The only difference between the two kinds of 
radiation is that the former is described by a spin-2 field, while the latter is 
spin one. As regards the two-level atom it is sufficient for us to assume that it 
has a non-zero cross section a for the absorption of gravitons at a frequency 
vq in its rest frame. The next point concerns the walls of the confining cavity. 
For the case of electromagnetic radiation it is safe to assume that one can 
construct walls with arbitrarily high absorptivity, and that the leakage time 
t L can be pushed to infinity. For the gravitational case, we assume that the 
cavity has an absorption coefficient f(v) = 1 — r)(v), so that for a spherical 
cavity of radius R 

d£{v)/dt _ r](u)p(u,t)cAnR 2 
E{y) ~ 47rR 3 p(u,t)/3 

where S(u) is the energy in the cavity at frequency v. 
This may be used to define a leakage time-scale 



(3) 



t L {y) w R/cqiy), (4) 
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and a frequency averaged scale 1l = R/cr], r] being the mean 'leakage' co- 
efficient. Clearly, any reference to the thickness or density of the absorbing 
material will be through r)(u). As expected, leakage takes more time for 
larger cavities. 

Consider next, the motion of the atom of mass M, when it is pushed about 
by the gravitational radiation, which has a time-dependent spectrum p(y,t). 
We assume that this motion can be described by the Langevin equation 
(2). This implies that the velocity of the atom will acquire a Maxwellian 
distribution over a time-scale M/a(z/ ), irrespective of its initial velocity. 
When this has happened for all possible frequencies, we can say that the 
radiation has also reached equilibrium. The time-scale t# should thus be 
defined as M/a m i n , where a m i n is the minimum value of a(v), obtained by 
varying v. A major task is to relate a min to an initial distribution po{y) and 
we address this question below. Right now, we compare i B with the average 
time-scale 1l defined earlier. Requiring ts <Jl gives 

R > Mcr]/a min (5) 

as the necessary condition for equilibrium. 

We must ensure that the box is not dense enough to undergo gravitational 
collapse. In other words, rj is assumed to be close to one. We also need to 
ensure that when R is chosen as in (5), the radiation in the cavity does not 
collapse to form a black hole. This requires 



2G.^i? 3 4 < Rc 2 (6) 
6 cr 



i.e., (ignoring numerical coefficients), 

R 2 < c 4 /Gp (7) 



p being the mean density of radiation. The compatibility of (5) and (7) 
requires 

1 /-2,.2 M 2 



M 2 r 2 n 2 

^-P- <R 2 < c^/Gp (8) 



or 



which may be taken to imply that p should satisfy 

P < * 2 mm c 2 /GM 2 V 2 . (9) 
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We interpret the equations (5)-(9) as follows. If the mean density satisfies 
the bound (9), then the radiation can be made to achieve equilibrium by 
choosing R within the bounds given by (8). The situation is different from 
the one for electromagnetism because r\ cannot be made one for electromag- 
netic radiation, and because we do not want t# < Jl to imply gravitational 
collapse. Such a constraint on rj is only to be expected, because unlike for the 
electromagnetic case, the principle of equivalence a priori fixes the 'charge to 
mass' ratio for gravity. 

It is evident that perfect absorption (/ = 1) is a sufficient but not a 
necessary condition for thermalization. This allows the construction of the 
above mechanism for achieving equilibrium. On the other hand, Smolin was 
mainly concerned with measurement of the quantum state of the gravitational 
field, and that certainly requires / ~ 1. Thus our results should be considered 
complimentary, rather than contradictory, to those of Smolin. 

We end this discussion by a comment on the coefficient a m i n in (9). Near 
equilibrium, a m i n {v) is proportional to is, and vanishes for v = 0, which only 
means that no momentum is absorbed when the energy levels coincide. So it 
is reasonable to define a m i n as a{u min ) where u min is the smallest frequency 
for which we are interested in studying the approach to equilibrium. Clearly, 
this issue arises for electromagnetic radiation also, and is not peculiar to 
gravity. 

For the equilibrium case, a{y) is simple to calculate; but for an arbitrary 
configuration an exact calculation is not possible. A very crude estimate 
for a m i n may be obtained as follows. Consider an atom which is moving 
with a speed V in an isotropic radiation with mean density ~p. The average 
momentum transferred to the atom per unit time will be of the order 

p V V 
-. — .co = pa. — 

c c c 

a being the mean cross-section of the atoms to absorption of gravitons 
gives a min ~ pcr/c, and from (9) it follows that 

p> GM 2 V 2 /a 2 . (11) 

Numerical estimates of the allowed ranges for R and p will depend on 
a, and on the nature of the confining box (through if). In principle, how- 
ever, gravitational radiation with p exceeding the bound in (11) will reach 
equilibrium if placed in a suitable box. 
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Having concluded the central part of our discussion, we would like to com- 
ment on several aspects of the analysis which deserve careful understanding. 

We begin with the key technique used in the analysis, viz. the Langevin 
equation and its validity. Langevin equation arises in many situations deal- 
ing with stochastic forces (see e.g. ,[6] for a discussion of its validity and 
comparison with Fokker-Planck equation) and is a good assumption if the 
scattering process is uncorrelated and homogeneous. In this particular case 
we assume that the system is interacting with a bath of gravitons which can 
be described classically as a stochastic background of gravitational waves. 
While one could think of more complicated scenarios, the description based 
on Langevin equation is adequate for our purpose. Of course, Langevin, 
equation is never "fundamental" and is more like a phenomenological equa- 
tion such as, say, Hooke's law. But since we are discussing a specific situation 
well within the domain of its validity, our results are trustworthy. 

The second, and somewhat more subtle issue is whether a depends on 
M and if so, how. This is interesting because, if a is proportional to M 
with, say a = Mao(u) where ao(u) is independent of M, then some of our 
bounds (which depend on a/M) will be independent of M. In fact, one may 
naively imagine that principle of equivalence will require a oc M making 
final results independent of M. In reality, however, a oc M 2 and there are 
different ways of seeing this. The first, and the simplest, is to note that 
in the electromagnetic case a oc q 2 (where q is the charge) and not a oc q; 
this is obvious since Langevin equation will treat both positive and negative 
charged particles at the same footing and the sign of a cannot depend on the 
sign of q. By analogy, the a for gravity should scale as (GM) 2 . The second 
argument is based on the result in equation (10) where we showed that a oc a. 
Now the lowest order scattering of gravitational waves by matter consists 
of matter oscillating under the infulence of the incident gravitational wave 
and radiating out gravitational waves due to its motion. This is completely 
analogous to the Thomson scattering of electromagnetic radiation and the 
scattering crosssection for this process can be worked out in a straight forward 
manner (Padmanabhan, 1978; unpublished). The result is as expected with 
a oc (GM/c 2 ) 2 ; i.e, the Schwarzschild radius acts as the analogue of classical 
electron radius. Once again, we find that a oc a oc M 2 . Finally, note that the 
averaging of the stochastic force acting on the particle will involve the two 
point correlation function of the force at the lowest order. That is why a will 
scale in proportion with the power spectrum [energy density] of the stochastic 
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background (quadratic in the gravitational wave amplitude) and — by the 
same token — will scale as M 2 . This imples that the time taken to reach 
equilibrium depends on the mass M of the atom. One should not use principle 
of equivalence naively to argue otherwise; principle of equivalence has no 
direct role to play in this issue. This is the same result as in the treatment of 
Brownian motion via the Langevin equation. As a result, the various bounds 
obtained above depend on M. Of course, the mass dependence of a is not 
germane to the main issue we are addressing. Our main concern here is to 
demonstrate that equilibrium can in principle be reached, even though it 
implies that the bounds will depend on M. 

It may appear that Eqns. (9) and (11) give conflicting bounds, the former 
being an upper bound and the latter a lower bound, on the density. What 
is happening here of course is that one does not in general know what a min 
is, and (11) results from using a very rough estimate for a min . It is also 
important to note that, we have tried to obtain fairly general conclusions 
with most plausible assumptions. In specific contexts, one may be able to 
strengthen the bounds or supercede them with other considerations but these 
will not be generic. For example, we have treated gravitational waves as a 
small perturbation on a flat background. While the gravity waves being a 
perturbation defines the context of the problem, one can think of the effects 
of background being curved. If the scale length of background curvature is 
comparable to the mean free path of gravitons, there will be further effects. 
We suspect that any such effect, which involves scattering of gravitons by 
the curvature will only go towards randomising the processes and will not 
affect our conclusions but this claim is hard to establish. 

Finally, we would like to speculate on some implication our results might 
have for the issue of information loss in black hole evaporation. It is well- 
known that a black hole radiates a Planck spectrum for linearized gravitons 
also, besides the Planck spectrum for other quantum fields of lower spin. 
For all practical purposes the radiated particles behave like a set of particles 
thermalised by some physical process. On the other hand, complete evapo- 
ration of a black hole resulting in purely thermal radiation as relic will lead 
to well-known information loss problems and possibly problems with unitary 
evolution of quantum mechanical systems. 

In the context of our work one could not help speculating whether the 
restrictions on thermalization of gravitational radiation by "normal" systems 
have any implications for the radiation of thermal gravitons by black holes. If 
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the black-hole graviton system can be treated like the 'gravitational radiation 
in a box' system discussed above, and if the system does not satisfy the con- 
ditions necessary for thermalization, the escaping non-thermal gravitational 
radiation could carry information. 

It is not clear whether one could interpret the black hole radiation as aris- 
ing out of a "thermalization" process; but if it is, then the possibility opens 
up that the gravitons can contain at least part of the information lost when 
the event horizon has formed and may lead to some interesting implications 
for the information loss problem. This idea, though rather speculative at this 
stage, deserves to be investigated further. 
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